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ABSTRACT 


A  representation  of  nonstationary  narrowband  random  processes  in  terms  of 
nonstationary  quadrature  components  is  proposed  in  a  form  analogous  to  that  used 
to  represent  wide  sense  stationary  narrowband  random  processes.  The  represen¬ 
tation  is  then  applied  to  a  specific  case  in  which  the  nonstationary  narrowband 
random  process  is  generated  by  the  product  of  white  noise  and  a  deterministic 
periodic  signal  and  .then  is  processed  by  a  narrowband  filter.  This  representation 
is  used  in  the  modeling  of  a  bi-level  pulsed  noise  jammer  which  is  assumed  to 
be  present  in  a  communication  channel.  The  effect  of  such  a  jammer  on  a  direct 
sequence,  binary  phase  shift  keyed  (DS-BPSK)  spread  spectrum  communication  re¬ 
ceiver  is  evaluated  and  characterized  in  terms  of  the  error  rate  performance  of  the 
receiver.  Families  of  performance  curves  are  plotted  to  demonstrate  the  effect  of 
various  parameters,  namely  signal-to-noise  ratio,  jammer  power  to  signal  power  ra¬ 
tio,  and  processing  gain,  on  the  error  rate  of  the  complete  spread  spectrum  receiver. 
The  analysis  carried  out  differentiates  between  two  cases,  namely  fast  jammers  and 
slow  jammers.  However,  the  analytical  tools  developed  make  it  possible  to  consider 
either  one  of  the  two  cases  without  resorting  to  quasi-stationary  arguments  as  has 
been  done  in  the  past. 
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I.  INTRODUCTION 

Electronic  Counter  Measure  (ECM)  techniques  have  gained  increasing  atten¬ 
tion  in  military  communications.  It  is  often  necessary  and  important  to  be  able  to 
analyze  and  predict  the  performance  of  communication  receivers,  especially  when 
operating  under  hostile  electronic  interference.  In  order  to  improve  the  survivabil¬ 
ity  of  tactical  communication  systems  operating  in  jamming  enviomments,  spread 
spectrum  modulation  techniques  have  been  introduced,  analyzed  and  implemented 
in  many  systems  in  order  to  mitigate  the  effects  of  intentional  electronic  inter¬ 
ference.  In  many  of  the  analyses  carried  out  to  date,  the  jammer  interference 
is  typically  assumed  to  be  a  stationary  (or  wide  sense  stationary)  random  process 

with  well-specified  statistical  characterizations  [Ref  1],  and  the  natural  interference 

* 

introduced  in  the  communication  channel  is  generally  modeled  as  additive  white 
Gaussian  noise  (AWGN).  However,  in  practice,  there  are  many  communication  en¬ 
vironments  which  do  not  fit  the  model  just  described.  Consider,  for  example,  a 
communication  system  which  is  being  jammed  by  a  continuous  wave  (CW)  tone 
near  the  transmitter  operating  center  frequency,  or  by  a  distorted  retransmission 
of  the  transmitter’s  own  signal.  The  interference  cannot  be  accurately  modeled 
as  a  stationary  random  process  in  either  case.  Another  typical  jamming  scenario 
involves  a  jamming  signal  that  may  be  pulsed  between  various  power  levels  at  a 
rapid  rate  and  is  therefore  not  stationary  [Ref  2]. 

The  effects  of  nonstationary  interference  on  the  performance  of  a  spread  spec¬ 
trum  communication  system  is  not  known  in  general.  It  is  however  possible  in 
many  cases  to  characterize  nonstationary  interferers  and  determine  their  effect  on 
specific  communication  systems.  In  fact,  if  the  statistics  of  the  nonstationary 
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process  are  periodic,  such  a  prc:ess  is  said  to  be  cyclostationary  [Ref.  3].  The 
statistics  and  power  spectrum  of  a  cyclostationary  process  are  normally  computed 
by  using  phase  randomization  or  time  averaging  techniques  [Ref.  3]. 

In  view  of  the  importance  of  the  use  of  nonconventional  forms  of  generating 
channel  interference,  the  aim  of  this  thesis  is  to  analyze  the  effects  of  nonstationary 
narrowband  interference  on  the  performance  of  a  spread  spectrum  communication 
receiver.  More  specifically,  the  results  of  the  analysis  are  applied  to  the  special  case 
of  a  narrowband  bi-level  pulsed  noise  jammer,  since  it  is  an  effective  form  of  ECM 
[Ref  3].  Although  certain  results  concerning  the  same  type  of  jammer  considered 
here  can  be  obtained  by  other  methods  [Reference  1  and  2]  the  results  derived  are 
based  on  the  application  of  quasi-stationary  assumptions  which  oftentimes  are  not 
realistic  and  cannot  always  be  justified  in  practice. 

In  Chapter  II  of  the  thesis,  methods  analogous  to  those  used  in  the  repre¬ 
sentation  of  narrowband  random  processes  are  used  to  characterize  a  stationary 
narrowband  nonstationary  random  process.  Such  a  random  process  is  shown  to 
be  specified  in  terms  of  its  nonstationary  quadrature  components.  The  expression 
of  the  autocorrelation  function  of  the  nonstationary  narrowband  random  process 
is  then  established,  and,  all  the  results  are  used  to  specify  the  characterization 
of  a  bi-level,  narrowband  pulsed  noise  jammer  which  can  be  obtained  by  ampli¬ 
tude  modulation  of  white  Gaussian  noise  which  is  filtered  by  a  narrowband  filter 
centered  at  some  high  frequency  /0 . 

In  Chapter  III,  analysis  was  carried  out  in  order  to  characterize  the  output 
of  the  spread  spectrum  receiver  consisting  of  a  despreader,  a  bandpass  filter  and 
a  correlation  detector,  when  a  nonstationary  narrowband  process  is  applied  at  the 
input  of  such  a  receiver.  This  characterization  is  obtained  in  terms  of  the  mean 


and  variance  of  the  Gaussian  random  variable  at  the  output  of  the  receiver’s  deci¬ 
sion  circuitry. 

Chapter  IV  provides  a  detailed  analysis  of  the  performance  of  a  spread  spec¬ 
trum  receiver.  The  binary  phase  shift  keyed  (BPSK)  modulation  scheme  was  se¬ 
lected  as  the  method  by  which  digital  information  is  transmitted  and  the  method 
of  spreading  the  signal  bandwidth  was  chosen  to  be  a  direct  sequence  spread  spec¬ 
trum  modulation  approach.  Noise  and  jamming  were  assumed  to  be  sources  of 
channel  interference,  using  the  models  developed  in  Chapter  3  for  their  character¬ 
ization.  The  receiver  performance  specification  was  accomplished  by  deriving  the 
error  probability  or  the  bit  error  rate.  The  error  performance  was  plotted  versus 
various  factors  that  take  into  account  signal  power,  interference  power  and  band¬ 
width  spreading.  The  effect  of  the  various  parameters  on  the  receiver’s  performance 
was  analyzed  and  discussed. 

From  the  general  representation  of  the  narrowband  nonstationary  random  pro¬ 
cess,  other  forms  of  nonstationary  interference  can  be  proposed  and  their  effect  on 
communication  receivers  can  be  derived  by  similar  approaches.  Although  the  de¬ 
rived  performance  of  the  spread  spectrm  receiver  was  developed  for  BPSK  signal 
modulation  with  direct  sequence  bandwidth  spreading,  the  general  results  can  be 
modified  and  extended  to  obtain  performance  evaluations  of  receivers  in  which 
other  modulation  schemes  are  employed. 


n.  NONSTATIONARY  RANDOM  PROCESSES 


A.  NONSTATIONARY  NARROWBAND  RANDOM  PROCESSES 

A  random  process  is  said  to  be  narrowband  if  its  power  spectral  density  is  zero 
except  for  a  narrow  frequency  region  (W  Hz)  around  a  high  carrier  frequency  (/„), 
that  is, 

W  «  2 nfa  and  2irf0  »  1 

If  X(t)  is  a  sample  function  of  such  a  real  random  process  with  zero  mean,  it  may 
be  expressed  as:  [Ref.  4] 

X(i)  =  Xe(t)  cos2irf0t  -  X,  (t)sm2irf0t  (2.1) 

where  the  power  spectral  density  of  the  random  processes  Xe(t)  and  X,  (t)  is  zero 
except  for  a  narrow  frequency  region  (W  Hz)  around  /  =  0.  The  random  processes 
Xe(t)  and  X.  (t)  are  referred  to  as  the  quadrature  components  of  X(t). 

Random  processes  can  be  generally  classified  into  stationary  and  nonstation¬ 
ary  processes.  If  the  joint  probability  density  function  (p.d.f.)  of  a  process  X{t) 
taken  at  arbitrary  times  tx,t7  ...t„  is  invariant  to  arbitrary  time  shifts  for  any 
integer  n,  such  a  stochastic  process  is  said  to  be  stationary  [Ref.  5].  A  random 
process  is  nonstationary  if  it  does  not  have  the  above  properties.  A  weaker  form 
of  stationarity  is  the  so-called  wide  sense  stationarity  (W.S.S.)  involving  first  and 
second  order  moments  of  a  random  process  [Ref.  6l.  The  representation  of  Eq.  2.1 
normally  implies  a  W.S.S.  random  process  X(t). 

We  propose  here  to  represent  a  nonstationary  narrowband  random  process  in 
a  form  analogous  to  that  of  Eq.  2.1  and  to  study  the  statistics  of  the  associated 
quadrature  components.  This  is  accomplished  by  analyzing  the  output  of  a 


narrowband  filter  when  the  input  is  an  arbitrary  nonstationary  process.  Before 
attempting  this,  a  convenient  method  of  specifying  a  narrowband  filter  is  developed 
first. 


1.  Lowpass  Equivalent  of  a  Narrowband  Filter 

Any  bandpass  filter  can  be  specified  in  terms  of  low-pass  filter  equivalents. 
Let  hL  (t)  be  the  impulse  response  of  a  low-pass  filter  and  HL  (/)  be  the  Fourier 
Transform  of  hL(t),  that  is, 

hL  (!)  <£>  Hl (/)  =  (2.2) 

where  HL  (/)  is  not  necessarily  symmetric  about  /  =  0  so  that  hL  ( t )  is  not  neces¬ 
sarily  real.  Then  the  Fourier  Transform  HB  (/)  of  a  bandpass  filter  having  impulse 
response  hB  (t)  can  be  specified  as  follows:  [Ref.  7] 

hB  (i)  Hb  (/)  =  HL  (f-  /.)  +  Hl  (-(/  +  /.))  (2.3) 

As  shown  in  Figure  2.1,  HB  (/)  is  realized  by  summing  the  frequency  up  shifted 
(by  fo)  and  the  frequency  reversed  and  down  shifted  (by  /„)  Fourier  Transform 
of  the  low-pass  filter.  The  negative  sign  in  the  term  HL  (— (/  +  /„))  is  due  to  the 
frequency  reversion  of  HL  (/)  described  above. 

We  can  remark  that  in  general  hL  (t)  is  a  complex  quantity.  However,  since 
the  focus  is  on  generating  real  narrowband  random  processes,  we  will  require  that 
hB  (t)  is  a  real  function  of  t.  By  inverse  Fourier  Transformation  of  Eq.  2.3  and  by 
expressing  hL  (t)  in  the  following  complex  form, 

hL(t)  =  hLr{t)  +  jhLi(t)  (2.4) 
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where  both  hLr(t)  end  hLi(t)  ere  real  functions  of  t,  we  can  obtain  the  impulse 
response  hB  (t)  of  a  narrowband  filter  as  [see  Appendix  A] 


hB  (0  =  2 Re  {hL  (0e,9'/** }  (2.5) 

where  i2e{(-)}  represents  the  real  part  of  the  function  in  brackets.  It  is  demon¬ 
strated  in  Appendix  A  that  hLr(t)  must  be  an  even  function  and  hLi(t)  must  be 
an  odd  function  if  hh(t)  is  to  be  real. 

2.  Output  of  a  Narrowband  Filter 

.In  the  case  of  a  stationary  random  process  input,  the  statistics  of  the 
output  of  a  narrowband  filter  is  well-known  [Ref.  4].  We  will  investigate  the  more 
general  case  where  the  input  random  process  to  a  narrowband  filter  is  real  but 
nonstationary.  Let  Y(t)  be  the  output  due  to  a  real  input  signal  X(t)  applied  to 
a  narrowband  filter  having  transfer  function  Ha  (/)  as  described  in  Section  1  [Fig. 
2.2].  Then  Y(t)  can  be  expressed  as  the  convolution  of  Jf(t)  and  hB  (f)  namely 

Y[t)  =  J~  hB  (t  -  a)X(a)da  (2.6) 

It  is  demonstrated  in  Appendix  B  that  Y (t)  can  be  further  expressed  as 

y(t)  =  2Xe(t)coa2irf0t  -  2X,(t)  sin2x/0t  (2.7) 


where 


Xc(t)=  r  [kLr(t 

*  —  oo 

*.(<)  =  r 

J  —  oo 


-  a) cos 2 xf0a  -I-  hLi(t  -  a)sin2x/0a]  X(a)da 

-  a)  cos 2 irf„a  -  hLr{t  -  a)sin2ir/0a]  X(a)da 


(2.8) 

(29) 
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There  is  a  certain  degree  of  similarity  in  the  mathematical  form  of  Eq.  2.7  and 
Eq.  2.1.  However,  the  properties  of  Xe[t)  and  X, (t)  given  by  Eq.  2.8  and  Eq.  2.9 
respectively  must  now  be  established. 

3.  Output  Correlation  Function  for  Nonstationarv  Input 

By  definition,  the  autocorrelation  function  Ry{ti,t3)  of  any  random  pro¬ 
cess  Y  (t)  is  the  expected  value  of  the  product  of  Y  (tx )  and  Y  (t3 )  for  arbitrary  time 
instants  fx  and  t3  [Ref.  5],  that  is 

R,(t1,ta)  =  ^{y(ii)y(^)}  (2.10) 


In  order  to  obtain  Ry  (tx ,  t2 )  for  the  output  Y  (t)  of  the  filter  described  in  Section 
2,  it  will  be  necessary  to  evaluate  the  following  quantities  (since  Y  (t)  is  the  linear 
combination  of  Xe(t)  and  X,  (£)) 

=  E{X.(h)X'(t7)}  .  (2.11) 

=  £{*.(«.)*.(«*)}  (212) 

=  E{X.[t,)X.[t,)}  (2.13) 

=  E{X.(h)X,(h)}  (2.14) 

The  four  quantities  above  are  complicated  functions  of  tx  and  t3  [See  App.  Bj.  It  can 
be  observed  that  in  general  Re{ti,t3)  #  and  Re,(ti,t2)  #  —  R,e(ti,t2); 

unlike  the  case  of  stationary  narrowband  process,  both  the  quadrature  compo¬ 
nents  have  identical  autocorrelation  function  and  the  two  cross-correlation  func¬ 
tions  are  null  if  both  Xc(t)  and  X,  ( t )  are  real  [Ref.  4].  The  autocorrelation  function 
Ry(tx,t2)  of  the  nonstationary  narrowband  process  Y(t)  can  be  formulated  as  fol¬ 
lows  [  See  App.  B]: 


Ry{t  li^)  —  ^ 


,  t2)  cos27r f0tx  cos  2ir  f0t2  +  R,  (t x ,  t2 )  sin 2nfa t  sin  nf0 12 


—  Re,  (tx,t2)  cos 2nf0tl  sin2irf0t2  —  R,e(tltt2)sin2irf0tcos27rf0t2 


(2.15) 
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It  is  clear  that  Ry(tl ,  t3 )  is  also  a  linear  combination  of  Re(t1,t2),  R,(tlit2),  Ret  {t  j ,  t2 ) 
and 


B.  A  MODEL  FOR  THE  NONSTATIONARY  INPUT 

The  results  of  the  previous  sections  will  be  applied  now  to  a  specific  general 
class  of  nonstationary  random  processes  which  are  generated  by  a  system  to  be 
described  now.  This  is  done  in  order  to  continue  the  development  for  a  class  of 
nonstationary  random  processes  of  interest  and  to  be  able  to  apply  the  results 
to  some  practical  problems  involving  jamming  of  spread  spectrum  communication 
systems.  Let  the  nonstationary  input  JY(t)  described  previously  be  the  result  of 
mixing  a  stationary  random  process  W  (t)  with  a  deterministic  signal  q(t)  as  shown 
in  Fig.  2.3.  We  further  assume  that  the  autocorrelation  function  Rw  ,  t2)  of  the 
stationary  random  process  W  (t)  is 

Rw(tl,t2)  =  E\r(tx)W{U)]  =S(t2 -t,)  (2.16) 


That  is,  W (t)  is  assumed  to  be  a  white  process  with  unit  power  spectral  density 
(PSD)  level.  Therefore 


R.(t1,t3)  =  iE(A'(t1)X(t3)j 

=  E[W(t1)?(tl)W(t3)g(f3)l 

=  03(fi)^(fa  —  *1)  =  72(ta)^(t3  —  ti) 


(2.17) 
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Pifura  2.3  Nonatationary  Input  X(t) 
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since  J2,  (tt ,  ta)  is  zero  unless  t2  =  tx .  We  can  now  apply  the  results  from  Appendix 
B  in  order  to  characterize  the  output  of  a  narrowband  filter  when  its  input  is  a 
nonstationary  random  process  of  the  type  being  described  in  this  section.  It  is 
demonstrated  in  Appendix  C  that  the  autocorrelation  terms  of  Equations  2.11  - 
2.14  simplify  significantly  when  X(t)  =  q(t)W(t),  namely 

(2-18) 

and 

(2.19) 

Finally,  the  autocorrelation  function  Ry[ti,t2)  for  the  nonstationary  process  F(t), 
namely  the  output  of  the  narrowband  filter,  has  the  following  form:  [see  App.  C] 

R*  (<1  >  *3 )  =  2  cos  2j r/0  (t3  -  ti )  f  hL  (tt  -  a)hL  (t3  -  a)q.  (a)da  (2.20) 

J  -  OO 

where 

«.(*)  —  **(*)  (2-21) 

Observe  that  Ry  (tx ,  t2 )  is  not  a  function  of  the  time  difference  (t3  —  ^ ),  so  that 
Y  [t)  is  (as  expected)  a  nonstationary  random  process. 

C.  A  SPECIFIC  EXAMPLE 

In  the  previous  sections,  we  have  established  a  representation  in  general  form 
for  a  narrowband  nonstationary  random  process  generated  by  applying  the  mixture 
process  q{t)W  (f)  at  the  input  of  a  narrowband  filter.  As  a  specific  example,  it  is  now 
desired  to  use  these  results  to  model  and  characterize  a  pulsed  noise  jammer  that 
is  processed  by  a  narrowband  filter.  Assuming  that  we  have  a  bi-level  pulsed  noise 
jammer,  we  can  model  the  jamming  signal  X{t)  as  the  result  of  mixing  white  noise 
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W(t)  of  unit  PSD  level  with  a  bi-level  deterministic  periodic  signal  q(t).  Such 
a  jamming  signal  present  in  the  communication  channel  would  cause  a  typical 
receiver  to  produce  a  jamming  component  at  the  output  of  the  receiving  front  end 
narrowband  filter.  This  component  is  denoted  Y (t)  as  shown  in  Figure  2.4. 

The  autocorrelation  function  R*  (fi ,  f3)  of  the  input  process  to  the  narrowband 
filter  is  given  by  Eq.  2.17  while  the  corresponding  output  autocorrelation  Ry  (fi ,  t2) 
is  given  by  Eq.  2.20.  T,  is  the  period  of  the  periodic  signal  q(t)  which  remains  at 
the  level  A  for  the  time  interval  [o,  pTq\  and  at  the  low  level  C  for  the  time  interval 
[pT,,T,], 0  <  p  <  1.  [Fig.  2.5]  The  narrowband  filter  HB  ( t )  is  based  on  a  low-pass 
equivalent  which  is  assumed  to  have  impulse  response  hL  (t)  =  B  •  e~BtU(t)  where 
{7(f)  is  the  unit  step  function,  and  B  specifies  the  3db  point  of  the  filter.  The 
details  of  the  derivation  are  presented  in  Appendix  D,  where  it  is  demonstrated 
that  for  the  special  case  under  consideration,  Eq.  2.20  becomes 

/^(f^f,)  =2cos27r/0(t3 
B(A3 

e2BT,  .e2BT,M  _e2BpT,  .  g3BTt[rm 

+  l-c-IBT,  J  / 

where  [ z ]  is  the  largest  integer  less  than  or  equal  to  z 


A- 


B  (<i  +  *i ) 


e3BUin(t„ 


**  >  ([r]  —  [«] 


(2.22) 


r  =  -p+  ^-Min(f1,ta), 


3  =  -1+  —  Min(tl,t3) 


(2.23) 

(2.24) 


and 


Min(tl,t3)  =  minimum  of  f,  and  f3 


Figure  2.4  Schematic  Representation  of  the  Bi-level  Jammer  Model 


Although  equation  2.22  is  not  a  simple  mathematical  form,  it  is  possible  to 
derive  the  average  power  of  the  random  process  Y(t),  by  time  averaging  Ry(t,t). 
The  details  are  again  presented  in  Appendix  D,  where  it  is  demonstrated  that 

JZ,(M)  =  B[A2p  +  c1{l-  p)\  (2.25) 

This  result  can  be  easily  verified  from  the  assumptions  made  at  the  beginning  of 
this  paragraph.  Furthermore,  it  can  be  seen  that  for  the  special  case  under  con¬ 
sideration,  Ry(t,t)  represents  the  jammer  power  at  the  output  of  the  narrowband 
filter  Hb  (/). 

It  has  thus  been  shown  that  a  nonstationary  narrowband  process  can  be  rep¬ 
resented  by  the  nonstationary  quadrature  components  in  a  form  analogous  to  that 
for  the  representation  of  wide  sense  stationary  narrowband  processes.  A  general 
expression  of  the  correlation  function  for  a  filtered  nonstationary  input  is  obtained. 
The  results  are  applied  to  a  special  case  in  order  to  model  a  bi-level  pulsed  noise 
jammer.  The  effects  on  a  spread  spectrum  communication  receiver  due  to  such  a 
pulsed  noise  jamming  signal  will  be  investigated  in  the  next  chapter. 


ITT.  EFFECT  OF  SPREAD  SPECTRUM  SYSTEM 
ON  NONSTATIQNARY  NOISE 

A.  MODEL  FOR  THE  SPREAD  SPECTRUM  RECEIVER 

In  spread  spectrum  modulation,  digital  information  is  transmitted  by  employ¬ 
ing  a  trasmission  bandwidth  which  is  much  larger  than  the  minimum  bandwidth 
required  to  transmit  the  digital  information  via  conventional  means.  In  the  spread 
spectrum  receiver,  the  demodulation  must  be  accomplished,  in  part,  by  correlation 
of  the  received  signal  with  a  replica  of  the  “spreading”  signal  used  in  the  transmit¬ 
ter  in  order  to  recover  the  digital  information  signal.  If  the  bandwidth  is  spread  by 
direct  modulation  of  a  data- modulated  carrier  with  a  wideband  spreading  signal 
or  code,  the  technique  is  referred  to  as  direct  sequence  (abbreviated  DS)  spread 
spectrum  modulation. 

The  simplest  form  of  DS  spread  spectrum  modulation  is  DS-BPSK  spread 
spectrum  transmission  where  a  binary  phase  shift  keyed  (BPSK)  modulated  car¬ 
rier  is  spread  by  a  wideband  “code  signal.”  [Ref.  2].  The  spreading  operation  can 
be  mathematically  represented  as  a  multiplication  of  the  carrier  by  a  function  c(f) 
which  takes  on  the  values  ±1  at  a  digital  rate  Re  periodically.  (Fig.  3.1)  The  trans¬ 
mitted  spread  spectrum  signal  is  received  together  with  some  type  of  interference 
and/or  Gaussian  noise.  Demodulation  is  accomplished  in  part  by  remodulating 
with  the  spreading  code  e(t).  The  process  is  commonly  referred  to  as  despreading 
and  is  a  critical  function  in  all  spread  spectrum  systems. 

A  simplified  model  of  the  DS-BPSK  receiver  is  illustrated  in  Fig.  3.2.  The 
figure  shows  the  received  signal  r(t)  first  processed  through  a  front-end  bandpass 


Figure  3.1  BPSK  Spreading  Code  c(t) 


Figure  3.2  Simplified  DS-BPSK  Spread  Spectrum  Reoeiver 


filter  centered  at  the  carrier  frequency  /0  and  having  a  bandwidth  wide  enough 
for  the  spreading  signal  (contaminated  by  thermal  noise  and  perhaps  a  jammer) 
to  be  allowed  through  with  negligible  distortion.  The  output  of  the  front-end 
bandpass  filter  is  then  mixed  with  a  local  synchronized  replica  of  the  code  spreading 
signal  c(t)  in  order  to  despread  the  signal.  The  despreading  operation  produces  a 
contaminated  BPSK  signal  from  which  the  digital  information  can  be  recovered  by 
conventional  techniques  using  a  filter  and  a  correlator  as  shown  in  Fig.  3.2. 

It  is  well-known  that  in  the  absence  of  any  jamming,  the  perfomance  of  the 
receiver  in  Fig.  3.2  in  terms  of  the  bit  error  probability  (Pt)  is  given  by:  [Ref.  2] 

*=«(vx)  (31) 

where  i2»  is  the  bit  energy  of  the  received  signal,  N„  is  the  one-sided  power  spec¬ 
tral  density  level  of  the  AWGN  interference  and  Q(-)  is  the  complementary  error 
function  defined  by 


B.  EFFECT  OF  NONSTATIONARY  NOISE  ON  THE  DESPREADER 

AND  BIT  DETECTOR  ON  NONSTATIONARY  NOISE 

The  jamming  signal  produced  by  the  bi-level  pulsed  noise  jammer  may  not 
necessarily  be  narrowband.  However,  after  the  filtering  operation  performed  by 
the  receiver’s  front-end  bandpass  filter  the  jamming  signal  can  be  considered  to 
be  narrowband.  We  therefore  investigate  here  the  effect  of  the  despreader  and  the 
bit  detector  system  on  the  input  nonstationary  narrowband  random  signal  which 
models  the  jammer. 

The  receiver  of  Figure  3.2  is  now  assumed  to  have  a  narrowband  nonstationary 
zero  mean  input  y(t)  as  illustrated  in  Figure  3.3.  The  output  of  the  bandpass  filter 
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Figure  3.3  Simplified  DS-BPSK  Spread  Spectrum  Receiver 
for  Nonstationary  Narrowband  Input 


is  denoted  by  *(t)  while  the  bandpass  filter  impulse  response  is  denoted  ht  (t)  and 
can  be  specified  in  terms  of  a  low-pass  filter  equivalents  denoted  hi  (t).  The  output 
of  the  receiver  due  to  the  input  y(t)  is  denoted  Yv .  It  can  be  seen  that 


*(<)  =  y(0«M  *  hi  (0  =  [  hi  (<  “  T)v(T)ciT)dr  (3-3) 

J  —  00 


and 

rT * 

Yj=  z(t)  cos  2irf0tdt 

Jo 

which  can  also  be  expressed  as 


-■/[/:  ht(t  -  r)y(T)c[T)di 


dt 


(3.4) 


The  jamming  of  the  transmitted  signal  is  always  being  considered  as  an  action 
of  adding  a  nonstationary  Gaussian  random  process  to  the  signal.  As  will  be  seen 
in  the  analysis  to  be  carried  out  in  Chapter  IV,  part  B,  the  effect  of  the  jammer  on 
the  receiver’s  performance  is  completely  accounted  for  by  evaluating  the  variance  of 
the  output  of  the  bit  detector.  That  is,  only  the  variance  of  Yj  will  be  established 
and  the  other  statistical  properties  of  Yj  which  have  no  direct  influence  on  the 
systems  bit-error  probability  will  be  neglected. 

The  variance  a]  can  be  shown  to  be  given  by  an  infinite  sum  of  terms,  namely 


Hi [f)Hi  {nRq  ~  f)Pe(-f)Pc(f  ~  nR ,)• 


f°°  Hl  (u)Hl  ( nR ,  -  u)  [5e(/  +  fo-u)  +  S.(f-f0-  «)]  dudf  (3.5) 

v-OO 

where  R .  is  the  period  of  the  nonstationary  jamming  signal  and  Cn  is  given  by: 


A3p  +  ca(  1  -  p) 
nrp 
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;  n  =  0 
J  n  £  0 


It  can  be  seen  from  Eq.  3.5  that  o  *  is  a  function  of  the  duration  Tb  of  the 
information  bit,  the  duration  Te  of  the  spreading  code,  chips,  and  Tf ,  the  jammer 
pulsing  period.  However,  practical  spread  spectrum  systems  utilize  spreading  codes 
for  which  the  chip  duration  is  very  much  smaller  than  the  information  data  bit 
duration  [Ref.  8].  That  is 


Te  «  Tb  or  Re  >>  Rt, 

where  Re  is  just  T~ 1 ,  the  chip  rate,  and  R+  is  Tk_  1 ,  the  bit  rate.  It  is  therefore 
possible  to  define  here  two  classes  of  jammers,  namely  fast  jammers  for  which  the 
prise  repetition  period  T,  is  of  the  same  order  of  magnitude  as  Te,  so  that  the 
jammer  pulsing  rate  is  much  larger  than  the  data  rate,  and  slow  jammers  for  which 
Te  is  of  the  same  order  of  magnitude  as  Tb.  The  effect  of  the  jammer  pulsing  rate 
on  the  variance  cr*  will  be  examined  separately  for  both  cases. 

1.  Variance  for  the  Fast  Jammer 

.Consider  first  the  Fourier  Transform  Pe(/)  of  the  function  Pc(/)  defined 
in  Eq.  E-2,  namely 

/>,(/)  =  /  cos2 =  G(/  -  /„)  +  G(f  +  /„)  (3.6) 

Jo 

where 

G(f)  =  J 

We  can  see  that  Pe{f)  is  significant  for  /  =  ±/0  over  a  restricted  frequency  band. 
Similarly,  Pe{(f  ~nR9)  is  significant  for  /  =  ±(/0  +nR,)  over  a  similarly  restricted 
frequency  band.  The  functions  Pe  (/)  and  Pc  (f  —  nRv )  are  diagramed  in  Figure  3.4 
from  which  it  is  possible  to  observe  that  if 

/o  nRq  +  Rb  <  <  fo  —  Rt, 
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or  equivalently,  if 


2  Rt  «  nR t 


(3.7) 


essentially  no  spectral  overlap  of  the  frequency  functions  occur.  This  condition 
except  for  n  =  0  is  satisfied  for  what  has  been  defined  as  a  fast  jammer  since 
R,  »  Rt  under  fast  jamming. 

It  is  possible  to  conclude  therefore  that  in  the  case  of  a  fast  pulsed  jamming 
environment,  the  only  contribution  due  to  the  infinite  sum  of  terms  that  make  up 
the  variance  a7  is  the  term  n  =  0  (See  Eq.  3.5).  Therefore 


*  =  c»  L 


r  HL(u)HL(-u)[S,(f  +  f„  —  «)  +  S.(/-/, -u))dudf 
J  -  oo 


(3.8) 


and  substituting  in  Eq.  3.8  the  corresponding  mathematical  description  of  the 
functions  Ht(-),HL(-),Pe( •)  and  S«(*)  as  given  in  Eqs.  E.3,  G.6,  G.8  and  G.17, 


results  in 

_  c0rk 

'  2 

8 

rM 

+  (2«- 

*r*  -  i) 

Be-‘ 

T,e-*r. 

,-BT, 


B9 


B{B9  +  *») 
bTk  [B2  -  6s)  bT,Te  [B7  -  b3  )7 


(3.9) 


v>n 


The  details  of  the  derivation  of  this  result  are  presented  in  Appendix  G. 

2.  Variance  o]  for  the  Slow  Jammer 

In  the  case  of  a  slow  pulsed  jamming  environment,  the  condition  R,  >> 
R «  is  not  satisfied,  therefore  a  certain  degree  of  overlapping  between  the  spectral 
factors  that  make  up  the  expression  for  a 9  will  occur.  While  the  infinite  sum  of 


Figure  3.4  Functions  Pc(f)  and  Pc(f-nRq) 


Eq.  3.5  can  be  approximated  accurately  by  a  finite  sum,  the  number  no ,  of  terms 
needed  in  the  approximation  will  be  such  that 


no-R,  >>  2R» 

or  equivalently 

«.  »  (3-io) 

In  this  case,  the  result  for  a *  becomes 

_  _  J—  oo 

r  ( u)Hl  ( kR ,  -  tt)  [Se  (/  +  f0  -  u)  +  S'  (/  -  fo  -  tt)]  dudf 

j  —  OO 

(3.11) 

It  is  also  clear  that  as  k  increases,  \Ck  |  decreases  so  that  the  sum  involved 
can  perhaps  be  further  truncated.  The  computation  of  a *  in  this  case  becomes 
quite  complicated,  however  using  numerical  integration,  a]  can  be  obtained  quite 
accurately. 


IV.  PERFORMANCE  OF  THE  SPREAD  SPECTRUM 
COMMUNICATION  SYSTEM 

A.  GENERAL  OUTLINE 

In  Chapter  III,  a  complete  derivation  for  o* ,  the  variance  of  the  Gaussian 
random  variable  that  quantifies  the  power  contribution  of  a  bi-level,  nonstationary 
narrowband  jamming  signal  at  the  output  of  a  DS-BPSK  Spread  Spectrum  Mod¬ 
ulation  receiver,  has  been  undertaken.  A  diagram  of  the  spread  spectrum  receiver 
is  shown  in  Fig.  3.1.  The  receiver  structure  is  optimized  based  on  principles  of 
Maximum  A  Posteriori  (MAP)  detection  for  binary  modulation  schemes.  The  re¬ 
ceived  signal  and  noise  is  first  processed  by  the  front-end  bandpass  filter  which  has 
a  bandwidth  just  wide  enough  to  pass  the  spread  signal  and  the  in-band  noise.  The 
narrowband  signal  is  then  processed  by  the  despreader  which  consists  of  a  multi¬ 
plier  mixing  the  received  signal  with  a  synchronozed  local  replica  of  the  spreading 
code  signal.  A  second  bandpass  filter  having  a  bandwidth  proportional  to  the  in¬ 
formation  bandwidth  follows  the  despreader,  so  as  to  further  eliminate  the  noise 
power  being  presented  to  the  remainder  of  the  receiver.  The  data  signal  is  then 
demodulated  by  a  correlator  and  finally,  a  decision  is  produced  every  Tb  seconds 
which  is  translated  into  0  or  1  data  bits.  This  generates  a  MAP  estimate  of  the 
digital  data  at  the  output  of  the  receiver. 

The  performance  of  such  a  receiver  is  usually  characterized  by  the  probability 
of  error,  which  is,  in  essence,  the  probability  of  an  incorrect  decision  made  at  the 
receiver.  The  probability  of  error  can  be  further  categorized  into  symbol  error  and 
bit  error  probability.  However,  for  the  case  of  Binary  Phase  Shift  Keyed  (BPSK) 


modulation,  the  two  error  probabilities  are  equivalent  and  provide  the  basis  for 
specifying  the  receiver’s  performance. 


In  the  next  section,  we  investigate  the  performance  for  the  case  in  which  a 
6PSK  modulated  signal  is  spread  by  direct  sequence  methods,  and  the  signal  is 
received  in  AWGN  interference  as  well  as  the  bi- level  pulsed  noise  jammer  described 
and  analyzed  in  previous  chapters. 


B.  RECEIVER  PERFORMANCE  FOR  BPSK  MODULATION 

The  structure  of  a  spread  spectrum  correlation  reciever  is  shown  once  again 

in  Fig.  4.1.  In  BPSK  modulation,  the  two  transmitted  signals  are  [Ref.  4] 

^2 ~Eb 

cos znj0i  u^e  ^  ib  ana 

(4.1) 


j,  (I)  = 


-o-V? 


cos2jt/0*  0  <  t  <  Tb 


where  Eb  is  the  energy  of  each  signal  and  Tb  is  the  duration  of  each  signal.  Since 
each  signal  represents  a  single  information  bit,  Eb  is  the  energy  per  bit  and  T~ 1  is 
the  bit  rate. 

At  the  output  of  the  second  bandpass  filter,  the  signal  y(t)  can  be  regarded  as 
the  summation  of  the  transmitted  signal  s<(t),  narrowband  noise  Ntb(t),  and  the 
jamming  signal  yj  (t)  convolved  with  the  impulse  response  ht  (t)  of  the  bandpass 
filter.  That  is, 


y(0  = 


*(0  +Nt„(t)  +yj(t)c(t) 


MO 


(4.2) 


*  =  1  or  2 


Hence, 

Y=  [  y(t)  cos27r  f0tdt 

Jo 


=  f  Si{t)  *  hj(t) cos 2irf0tdt  +  f  Nih(t)C(t)  *  h/(t) cos 2nf0tdt  (4.3) 
Jo  Jo 


or  simply, 


where  Sit  Ntk  and  Yj  represents  the  three  integrals  in  Equation  4.3  respectively. 
It  is  clear  that  Y  is  in  fact  a  random  variable  whose  statistics  can  be  evaluated 
once  certain  properties  of  5/ ,  Nth  and  Yj  are  established. 

We  first  assume  that  the  noise  contribution  at  the  front  end  of  the  receiver  is 
zero  mean,  white  Gaussian  noise  with  ^  as  its  two  sided  power  spectral  density 
level  [Fig.  4.2].  The  autocorrelation  function  for  the  noise,  ednoted  RN  (t,r)  is 

(4.5) 


Y  =  S,+N,h+  Y, 


(4.4) 


R«(t,r)=^S(t-T) 


where  6(f)  is  the  Dirac  Delta  function.  Based  on  these  assumptions  it  is  shown  in 
Appendix  H  that  Nth  is  also  a  zero  mean,  Gaussian  random  variable  with  variance 
o*h  specified  in  the  sequel.  Due  to  the  fact  that  the  receiver 'is  linear,  the  noise 
contribution  at  the  output  retains  its  Gaussian  statistics.  The  expression  for  a*h 
is  from  App.  H  given  by  (Eq.  H.18). 


°th  = 


NaTbk 
2x 


(4.6) 


k  =  0.903 


Re  ^  .  • 

Ge  —  —  =  processing  gam. 
Rh 


It  is  shown  that  and  S2  are  given  by  [Ref.  2] 

Si  =  J^r  and  s»  =  / 


EbTb 


(4.7) 
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Furthermore,  Yj  is  the  quantity  related  to  the  jammer  and  its  variance  a*  is  given 
by  (see  Eq.  H.4) 


CaTb 

8 


J(G.) 


where  J{Ge)  is  defined  by  Eq.  4.6.  To  further  simplify  the  notation,  let  X  = 
Nth  +  Yj  ,  so  that 

(4.7) 


_2  _  _2  ,  _3 
°X  ~  °th  +  ° J 


since  the  noise  and  jamming  are  zero  mean  uncorrelated  processes,  and 


Y  =  S{  +  X 


(4.8) 


then  X  is  zero  mean  Gaussian  random  variable  with  pdf  given  by  [Ref.  2] 


/.  (X)  = 


V-2 


:€  X 


not 


(4.9) 


For  the  BPSK  modulation  scheme  that  is  being  considered  here,  the  decision 
regions  of  the  signal  space  diagram  can  be  illustrated  as  shown  in  Fig.  4.3  [Ref  4[: 
From  equations  4.8  and  4.9,  it  can  be  seen  that  the  probability  of  receiving  the 
signal  Y  given  that  a  signal  S<  was  actually  transmitted  is: 


Pr{r  =  YlSi{t)  was  transmitted}  =  fr/,.  ( Y /Si )  =  fx(Y  —  Si) 


The  probability  of  bit  error  rate  is  then  given  by:  [Ref.  2] 


(4.10) 


P»  =  Q 


1 EbTb 

2ax 


or 


Q{\l  2«V<7?)) 


(4.11) 
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Figure  4.3  Signal  Space  Diagram  for  BKSP  Modulation 
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where  Q(-)  is  the  complementary  error  function  defined  earlier  in  Chapter  III. 
Substituting  and  a2j  from  Equations  H.4  and  H.18  results  in 


or 


(4.12) 


C.  ANALYSIS  OF  ERROR  PERFORMANCE 

As  derived  in  the  previous  section,  Eq.  4.16  yields  the  probability  of  bit  error 
Pb  for  a  DS-BPSK  spread  spectrum  receiver  under  AWGN  noise  interference  and 
narrowband  bi-level  fast  pulsed  noise  jamming.  Clearly  Ph  is  a  function  of  the 
power  spectral  density  level  N„  of  the  AWGN,  the  average  power  of  the  jammer 
given  by  BCa  .and  also,  the  processing  gain  Gc .  It  is  clear  from  Eq.  4.11  that 
the  receiver  bit  error  probability  will  increase  if  the  noise  and/or  jammer  power 
increases  as  the  Q(-)  function  is  monotonically  decreasing. 

In  order  to  analyze  the  quantity  Pb  in  more  detail,  we  define  the  signal  to  noise 
ratio  (SNR)  to  be  the  ratio  of  Eb  to  Na  and  the  jammer  power  to  signal  power  ratio 
(JSR)  to  be  the  ratio  of  average  jammer  power  Pj  to  the  average  signal  power  P 
[Ref.  4j.  Hence, 

SNR  =  ~  ,JSR  =  ?±  (4.17) 

2  V0  2 

where 


Pj  —  BCa  —  2irReC0  (4.18) 

P  =  EtR*  (4.19) 
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so  that  Eq.  4.16  can  be  reexpressed  as 

a  =  o( 

Based  on  Eq.  4.20,  families  of  curves  of  Pb  versus  SNR  for  fixed  JSR  values  can  be 
plotted  and  analyzed.  However,  first  plots  of  Pk  versus  JSR  for  different  values  of 
Ge  under  various  SNR  conditions  are  shown.  Logarithm  scales  are  used  for  Ph  due 
to  the  large  variation  in  this  parameter.  From  Figure  4.4a  to  d,  one  can  observe 
that  in  general,  under  intense  jamming  (i.e.,  JSR  large),  the  receiver  probability  of 
error  is  also  large,  thus  implying  that  a  significant  number  of  delivered  data  bits 
are  in  error.  Under  a  fixed  JSR  and  fixed  Ge  values,  it  can  be  seen  that  if  the  SNR 
is  large,  Pb  can  be  made  significantly  small.  For  example  in  Fig.  4.4b,  it  is  shown 
that  when  JSR  =  15dB  and  Ge  =100; 

Pk(at  SNR  =  4)  «  10"  3  7 

and 

Pk( at  SNR  =  20)  »  10- 30 

However,  when  the  jamming  is  very  strong  (for  JSR  >  30  dB)  Pb  tends  to  a  limit  of 
approximately  10“ 1  regardless  of  the  value  of  all  other  factors.  It  can  be  observed 
from  Eq.  4.16  that  if  JSR  increases  indefinitely,  the  limit  of  Pb  equal  to  0.5  is 
reached. 

From  Fig.  4.4a  to  d,  we  can  also  conclude  that  for  any  fixed  values  of  JSR 
and  SNR,  an  increasing  processing  gain  Gc  implies  a  decreasing  value  of  P„.  The 
phenomenon  is  more  prominent  in  a  high  SNR  environment  than  in  a  low  SNR 
environment.  There  is  no  significant  improvement  in  Pb  for  SNR  in  the  range  from 
0  dB  to  8  dB  however  a  large  improvement  in  Pb  is  obtained  when  SNR  exceeds 
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Figure  4.4  Pb  versus  JSR  for  (a)  Gc-10 
(b)  Gc= 100 
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Figure  4.4  Pb  versus  JSR  for  (c)  Gc-200 
(d)  Gc-400 
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12  dB.  The  immediate  conclusion  that  we  can  draw  from  this  is  that  by  increasing 
the  chip  rate  ( R ,)  of  the  spreading  code,  resulting  in  a  higher  processing  gain, 
the  perfromance  of  the  receiver  can  be  improved  under  the  jamming  environment 
analysed,  however  this  counter  measure  loses  effectiveness  when  powerful  jamming 
is  present. 

To  better  illustrate  the  effect  of  Ge  on  Pk,  another  family  of  curves  is  plotted. 
In  Fig.  4.5a  to  d  the  variation  of  Pk  as  a  function  of  JSR  is  shown  for  different 
values  of  G„  and  SNR.  It  is  clearly  shown  that  all  the  curves  converge  toward  Pb 
equal  to  0.5  as  JSR  increases.  From  Fig.  4.5a  we  can  see  that  Pk  (at  JSR  =  20 
dB)  reduces  from  approximately  10"  07  to  10" 14  when  Ge  increases  from  2  to  400 
under  weak  signal  (SNR  =  4  dB)  conditions,  while  Ph  reduces  from  approximately 
10" 0  7  to  10" 33  under  the  same  conditions  when  SNR  =  20  dB  (see  Fig.  4.5d). 
At  JSR  >  35  dB,  no  significant  gain  can  be  obtained  regardless  of  signal  strength. 
As  a  result  of  this,  we  can  conclude  that  by  increasing  the  chip  rate  to  counter  the 
narrowband  bi-level  fast  pulsed  noise  jammer  is  not  effective  under  strong  jamming 
conditions,  however,  a  high  chip  rate  can  reduce  the  probability  of  bit  error  of  the 
jamming  is  not  too  severe. 

It  is  also  possible  to  represent  the  error  performance  curves  by  plotting  log10  P* 
versus  SJR,  the  signal  power  to  jammer  power  ratio  under  various  SNR  conditions. 
Fig.  4.6  a  to  d  show  how  P»  varies  with  SJR.  The  curves  can  be  interpreted 
in  essentially  the  same  way  as  the  previously  mentioned  graphs.  By  comparing 
Fig.  4.6  a  to  d  with  the  results  obtained  in  Reference  2,  we  observe  that  tho 
performance  of  the  DS-BPSK  spreak  spectrum  reciever  degrades  under  a  jamming 
environment  in  a  manner  similar  to  the  case  of  barrage  noise  jamming  or  partial 
band  noise  jamming  environment  presented  in  the  reference.  Furthermore,  the 
receiver  performance  is  well  below  the  boundary  of  the  performance  for  the 
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Figure  4.5  Pb  versus  JSR  for  (•)  SNR  =  0  dB 
(b)  SNR*4  dB 


l  •  r ••  .r 


WO 


cv^v*: 


KRSSSKSKSi 


worst  case  pulsed  noise  jamming  derived  in  Reference  2  based  on  a  quasi-stationary 
analysis  method.  Except  for  the  weak  signal  condition  case,  the  receiver  perfor¬ 
mance  approaches  the  above  mentioned  boundary.  In  general,  by  increasing  the 
processing  gain,  the  performance  of  the  receiver  improves  significantly. 


V.  CONCLUSIONS 


In  this  thesis  it  has  been  demonstrated  that  any  nonstationary  random  process 
that  is  bandpass  filtered  can  be  represented  in  a  mathematical  form  analogous  to 
that  used  to  represent  narrowband  wide  sense  stationary  random  processes.  The 
quadrature  components  of  the  narrowband  nonstationary  process  are  themselves 
nonstationary  low  pass  random  processes.  In  general,  the  autocorrelation  function 
of  the  two  quadrature  components  are  not  identical  and  the  cross-correlation  func¬ 
tion  of  the  quadrature  components  is  not  necessarily  an  odd  function.  However, 
when  the  input  nonstationary  random  process  is  real,  the  correlation  functions  of 
the  quadrature  components  of  the  narrowband  output  process  have  properties  sim¬ 
ilar  to  those  associated  with  narrowband  wide  sense  stationary  random  processes 
The  general  representation  of  the  nonstationary  narrowband  random  process 
has  been  applied  as  a  model  for  a  pecific  class  of  electronic  counter  measure 
signals.  Specifically  a  bi-level  pulsed  noise  jammer,  in  which  the  power  of  the 
jammer  periodically  pulses  between  two  levels  is  modeled  using  the  previously 
mentioned  nonstationary  narrowband  random  process  representation  and  its  effect 
on  a  DS-BPSK  receiver  analyzed.  This  type  of  jammer  was  analyzed  because  it  is 
quite  effective  in  the  sense  that  it  produces  significant  performance  degradations 
on  the  receiver.  In  practical  cases,  the  jammer  pulses  between  two  levels  at  a  rate 
faster  than  that  of  the  data  bit  rate.  The  quasi-stationary  analysis  method  for 
receiver  performance  is  not  valid.  Therefore,  a  methodology  has  been  developed 
that  allows  rather  precise  evaluation  of  receiver  performance  in  terms  of  signal, 
noise,  and  jamming  powers  as  well  as  spread  spectrum  processing  gain. 
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The  effect  of  the  pulsed  noise  jammer  on  the  spread  spectrum  binary  phase 
shift  keyed  receiver  was  examined  in  detail  and  it  is  shown  that  a  strong  transmitted 
signal  is  more  effective  at  mitigating  the  effect  of  the  jammer  than  increasing  the 
chip  rate  of  the  spreading  code.  The  pulsed  noise  jamming  earn  be  shown  to  be 
quite  effective  under  various  signal  power  and  processing  gain  values,  as  rather 
large  bit  error  rates  result  when  powerful  jamming  is  present. 

The  thesis  develops  a  methodology  that  can  be  used  to  evaluate  the  effect  of 
a  nonstationary  narrowband  random  process  interference  on  communication  chan¬ 
nels.  There  are  two  benefits  from  this  study.  First,  in  terms  of  ECM  performance 
it  makes  it  possible  to  evaluate  the  effectiveness  of  such  jammers,  and  second,  in 
terms  of  electronic  counter  counter  measures  (ECCM),  the  performance  of  jam  re¬ 
sistant  receivers  can  be  evaluated  more  accurately.  Although  the  example  chosen 
in  the  study  involves  BPSK  modulation  spread  by  direct  sequence  methods,  results 
for  other  forms  of  nonstationary  jammings  on  different  types  of  spread  spectrum 
receivers  can  also  be  obtained  from  the  above  results  with  appropriate  modification. 
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APPENDIX  A 


REPRESENTATION  OF  NARROWBAND  FILTERS 


Let 


hL(t)^  HL(i)  =  A(iy,,n 
be  the  specification  of  a  low-pass  fitter  (LPF),  and  let 

As  (t)  ~  HB  (/)  =  Hl  (/  -  /„)  +  Hl  (-(/  +  /„)) 


(-4.1) 


(-4.2) 


where  /„  is  assumed  to  be  much  higher  than  the  cutoff  frequency  of  the  LPF. 
Thus,  Hb  (/)  is  the  transfer  function  of  a  narrowband  filter  centered  at  f  =  ±f0. 
By  inverse  Fourier  Transformation,  we  have: 

As(t)=  f  Hb  (f)e’2wItdf 
J-OO 

=  f°"  am/ -/>)'”•'’<*/+  /”  (A3) 

j  —  CO  J  —  oo 

=  hL(t)e’2n/-t  +hL{-t)ri2nf't 

If  we  express 

MO  =  hLr{t)  +  jhLi{t) 
since  hL  (t)  may  be  complex,  then 

hB  (0  =  [hLr  (0  +  M  (~01  cos  Z^fot  +  [hLi  (-t)  -  hLi  (0]  sin  27r/0t+ 

i  {[M(0  +  hLi(-01cos27r/ot  +  [M(~0  -  HLr(-t)\sm2xf0t}  (A.4) 

Since  hB  (t)  will  be  assumed  to  be  real,  the  imaginary  part  in  Eq.  A.4  must  vanish, 


so  that 


^Li  (0  —  (  —  0 


hLi  is  an  odd  function. 


(A.5) 


and 


hLr(t)  =  ^Lr{— t)  =>  hLr  is  an  even  function.  (.<4.6) 


Observe  also  that 


hL  (-0  =hLr(-t)+  jhLi  (-0  =  hLr  (t)  -  jhLi  (0  =  hi  (<)  (A.7) 


Therefore 


and  finally 


hB  (f)  =  2hir{t)  cos2ir/0i  -  2hZ(i(t)  sin2jr/0f 


MO 


(X.8) 


M.9) 
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OUTPUT  OF  A  NARROWBAND  FILTER 


The  output  y(t)  of  a  narrowband  filter  due  to  an  input  x(t)  is 

y(t)  =  J  hB(t-  a)x{a)da 
J  —  00 

=  f  [2 hLr  ( t  -  a)  cos  2 irf0 ( t  -  a)  -  2 hLi  (t  -  a)  sin  2 xf0 (t  -  a)]  x(a)da 

J  —  oo 

(B.l) 

where  the  representation  derived  in  App.  A  has  been  used  to  obtain  Eq.  B.l.  By 
expanding  cos  2x  fa  (t  —  a)  and  sin  2 n  fa  (t  —  a)  it  is  possible  to  obtain 

y(t)  —  2Xe(t)cos2Trfat  —  2X,  (t)  sin27r/0f  (5.2) 

where 

Xc(t)  =  j  [hLr{t  -  a)cos2irf0a  +  hLi(t  -  a)  sin27r/0a] X(a)da  (5. 3) 

J  -  OO 

and 

X,(t)  =  j  [hLi(t  -  a)coa2irf0a  -  hLr(t  -  a)sin2irf0a]X(a)da  (5.4) 

so  that  a  standard  narrowband-type  representation  is  possible  for  the  output  y(f). 
If  .X'(t)  is  a  random  process  that  is  not  necessarily  stationary,  Xe(t)  and  X, (f)  are 
themselves  random  processes  with  autocorrelation  functions 


=  j”  r  -<*)««"•• +fcl(«. 

v  “  OO  *  — •  OO 

•  \  R.(a,0)dadff 

4  f  oo  p  oo 

=  2  /  /  -/3)S2”'°la  +  <’)}R,(«,0)dad3+ 

J  —  OO  J  —  oo 


\  r  f°  Re  {/.I  (t, ,  a)hL  (t,  -  ny” }  R,  (aj)iad0  (B.  5) 

and 

*.(<!,<»)  =  E{X.(*,)X.(t,)} 

=  -\j”  f~  Re  {M<>  -0)'-‘”'-l-t,)}R,{<*,W<*dl3+ 

\T  f  Re{AL((,  -o)At(ls  -/3)e',"-|-'l}ft,(a,/S)<Ja<i(3 

*  “  OO  w  —  OO 

(B.6) 

The  cross  correlation  functions  are  given  by 
=E{X.((,)X.((,)} 

=  j  jH  J~  Im[kL{ti  -  a)hL{t,  -  R,(a,P)iadP+ 

\r  f°°  ImU-L(tl-a)hL(t,-0)ri,’M-*»\R.(a,0)dad0  (B. 7) 

*  V-OO  J  —  oo 

and 

E..(<1,iJ)  =  E{X.((I)X«(tJ)} 

4  /*  OO  /*  oo 

=  -  /  /  /m{h*  (fx  -a)M*a  -/?)e-^-(-'»)}JR,(a,/?)dadj9+ 

•  —  OO  *  *  OO 

-i/”  /”  /«{*!(*, 

*  —  oo  *  —  oo 

(B.8) 

Observe  that  in  general  Xc  (f)  and  X,  ( t )  normally  labeled  the  quadrature  compo¬ 
nents  of  the  narrowband  process  y(t)  are  not  stationary  in  general.  Nevertheless, 
the  derived  autocorrelations  and  cross  correlations  can  be  used  to  obtain  the  au¬ 
tocorrelation  of  the  random  process  of  y(f),  namely 
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.R»(*i.*a)  =  £{y(*i)y(fa)} 

=  E^2[Xe{tl)coe2nf0tl  -  X,{tl)sm2vf0tl] 

•  2[ Xe  (t7 )  cos  2xf0t2  -  Xt  (*3 )  *in  2irf„  t3  ]  | 

=  4  Jic(tlft2)  cos2irf0t1  cos27r/„f3  +  R.  (ti  ,t7)  sin27r/0t1  sin2jr/0ta-f 

-  ^.(ti  ,f3)  cos2ar/0ti  sin2jr/<,i3 

-  ^.«(<i»<a)8in2s-/0i1  cos2*/0f3]  (£.9) 


NARROWBAND  PROCESS  AUTOCORRELATION  FUNCTION  EXAMPLE 


The  autocorrelation  functions  determined  in  Appendix  B  are  now  used  to 
treat  a  specific  case,  namely  when  the  input  process  X(t)  to  the  narrowband  filter 
is  generated  by  the  system  of  Fig.  2.3.  The  resulting  output  y{t)  is  a  narrowband 
process  bearing  quadrature  components  whose  autocorrelation  function  is  given  by 
(see  Eq.  2.17) 


Rc(ti,t7)  =  ^  J  R e{hc(tl  -a)hL{t7  -  a)e+’2"  2/”“  }  q,  (a)da+ 

J  —  oo 

^  J  -  a)}  q.  (a) da 


(C.l) 


Given  that  hL  (•)  represents  a  low-pass  system  and  q,  (a)  would  normally  vary  much 
more  slowly  than  the  frequency  2 f„,  it  can  be  seen  that  integrals  involving  a  2 f0 
frequency  are  negligible  resulting  in  further  simplifications.  Thus 

Re(ti,t3)  =  j  J  Re{/i2(*i  -a)hL[t7  -a)}q.[a)da  (C. 2) 

*  —  OP 


and  similarly 


R*(ti  't*)  =  \f  R*  “  a)M*a  ~  a)}q,{a)da  =  Rc[tlft7)  (C.3) 

v  —  OO 


The  crosscorrelations  are  given  by 


Rc.{ti,t2)  =  \  f  Im{hL{tl  -  a)hL(t7  -  a)}q,{a)da  =  ( C.4 ) 


For  the  special  case  where  is  real,  J2..(ti,ta)  and  R,e{ti,t3)  are  both  equal 
to  zero,  so  that  from  Eq.  B.9  and  Eqs.  C.2  and  C.3,  Ry{t1  ,t3)  becomes 


R* (<i . <a)  =  4/2. (<i , t3 )  (cos  2rf0tl  cos  2 *f0t3  +  sin 2xf„ti  sin 2 wfat3 


=  4J2,  (tx ,  t3 )  cos  2 wf0  (t3  ~  ) 


(C.5) 


and  finally 


i2,(ti,ta)  =  2coe2jr/0(ta  -  tx 


(C.6) 
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AUTOCORRELATION  FUNCTION  AND  AVERAGE  POWER  AT  THE  OUT¬ 
PUT  OF  A  NARROWBAND  FILTER  DUE  TO  A  BI-LEVEL  PULSED  JAMMER 


When  ?(t)  is  a  bi-level  signal  as  shown  in  Fig.  2.5,  it  is  clear  that 

«.(<)  =  A1  —  j(A’  -  <’)  Y.  I1'!'  -  *7.  -  pT,)  -  £7(1  -  kT,  -  T,)|  J  (D.l) 

where  U(t)  is  the  unit  step  function.  Let 

I(t)  =  f  M*i  ~  a)^z.(«a  “  a)q.(a)da 

This  integral  must  be  determined  in  order  to  evaluate  as  given  by  Eq. 

C.6.  Substituting  q.  (t)  into  I(t)  and  utilising  the  one  pole  filter  specification  for 
hL  (t)  having  3  dB  cutoff  B  rpa,  we  get 


/(t)  =  +  _  B2(A9  -  C3)e~0,,|  +  '*) 

[|"e,1,'i;(fl-a)^(ta-o)  £  U(a  -  hT,  -  pT,)da 

fc  ■  “  00 

-  r  €>BmU{tl-a)U{t9-a)  Yl  U(a-kT,  -T,)da] 

W  —  00 


(D.  2) 


fc  *  —  GO 


The  remaining  integrals  in  Eq.  D.2  are: 


«  <*M  in  ( 1 1 .  t ,  )  .  l*i 

>  /  e3Bada=—  T  _  Mi»tnr, 


if  (fc  +  l)Tf  <  Min(ti  ,t3) 


-  0  otherwise 


(D.  3) 


and 


»  /.UiDftx.t,)  !  ill  r  1 

£  /  e3Bada=  —  £  [eaflM“(t*'<*)  _*»*<*+ or,  1 

hm-oo  •'(k_1)r«  °  Jks-oo 

if  (fc  +  l)T,  <Min(i1,«a) 

=  0  otherwise 

where  [z\  =  largest  integer  less  than  or  equal  to  z  and 


(0.4) 


r  =  -t  +  — Mm((,,la) 


S  =  1  +  ^-Min(tltt3) 
■*« 


(0.5) 


The  difference  of  the  two  integrals  given  by  Eq.  D.3  and  D.4  produces 


i  (  M  r 

55  E  t 
1*=-°° 


eJBMin(*1,»,)  _e2fl(*  +  o)T, 


i  l#l  r 

]-  E. 


^2B  M  in  (t  i  _  £B  ( fc+  1  )Tq 


km  —  oo 


(0.6) 


since  0  <  p  <  1,  we  have  r  >  s,  so  that 


^  e2BUiu(t,.t$)  _  £  [[r]  _[5]]  (£>.7) 


—  oo 


—  oo 


Furthermore 

H  jli,  -3«rt  -2BT,[.|  _  -3B4>rf  -2flTf[r| 

J2  e3B{k+l)T<  -  J2  e3B(k+')T'=- - - - - _  e_3-Tt - - -  (P.8) 


km  —  OO 


Combining  the  above  results  yields  (see  Eq.  C.6) 

X,(tltt3)  =  2 cos 2irf0t  •  /(*) 
n _ A  A?  & 


=  2coe2*/0t<  +  — 


—A~2~ lh-Bltl+,i)  eaflM‘“<‘l  *»)  ^[r]  -  [a]^ 

e2flr,e2flT,|.|  _  e3flpTf  g2BTf|r|  -I  1 

+  l-e-3Br.  Jj 


(P.9) 


1  -e~2BT' 


r.  ^  v*  r  i 


a 


<  i 

i 

*» 


I 

<7 


Since  X(t)  represents  a  bi-level  pulsed  jammer,  the  power  due  to  the  jammer 
at  the  output  of  the  narrowband  filter  denoted  by  Pj  can  be  obtained  by  letting 
*i  =  =  t  in  the  expression  of  Ry(ti,t7)  and  time  averaging  the  result.  Thus 


*,(<■<)  =  A'B  -  B(A7  -  C*)([r]  -  [*])  - 

* 

_  g2 BpT,  e-3B(«-T,|r|) 


-as  (t-r,|«|) 


_  e3BpTte-3B(«-r,[r|)j 

(£>.10) 

which  must  now  be  time  averaged  in  order  to  obtain  Pj ,  the  average  jammer  power. 
With  the  aid  of  Diagram  D.l,  D.2,  and  D.3  we  note  that  all  the  three  functions 
{[r]  —  [b]}  ,  {*  —  T,  [«]}  and  {t  —  Tf  [r]}  are  functions  of  t,  and  their  time  average 
values  can  be  calculated  as  follows: 


(M  -  [SD  =  Jim  (^r(r.  +  (n -  1)(1  - 
=  lim  -fp+  (1  -./>)n]  =  1  -p 

nr-  oo  Tl 


g-3B(t-Tf[S|) 


=  lim  f 
n-oo  nT„  J0 


1  /*r«  --2BT, 

4-  /  ne~2B{t+t,)  = _ 


2  BT, 


(£>.11) 


(1  —  c~2BT')  (£>.  12) 


,-3B(l- 


-  1  rPT'  fr* 

*-r«H>  =  lim  7 - /  c~2Btdt  +  n  iT2B<‘+ 

(n  +  ^r,  J0 


pT')dt 


e-3B,r,(i  _g-3BTt) 

“  2BT' 


(£>.13) 


Therefore,  the  final  expression  for  P}  is 


Pj  =  A2 B  -  B[A 2  -  C2)(l  -p)~ 
=  B[A2 p  +  C2(l  -  />)] 


B{A2  -  C2)  [1  -c~2BT'  l-e"2BT« 


X_e-3BTf  j.  2BT, 


2BT„ 


(*>•14) 


V  V  •  V  V  V  »*. 


V  *  .  •  i 


3 

2 


1 


-  — »  I 

(l-MTj)  (2-MTj) 


* 


f  W“[s] 


r  tTq 


&+*r») 


< 


Figure  D.l  Function  [r]-[s] 
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Function  (t-Tq[s]) 
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lPPENDDC 


VARIANCE  OF  Y, 


Since  y(t)  is  Gaussian  and  zero  mean  random  process,  Yj  is  a  Gaussian,  zero 
mean  random  variable  due  to  the  fact  that  linear  operations  only  axe  applied  to 
y(t).  [Ref  2]  The  variance  of  Yj,  denoted  by  a]  can  therefore  be  expressed  as 


follows: 


''"Iff  z(t)z(v)  cos  It  f0 1  cos  2t  f0  vdtdv| 


rT*  rTb 

—  /  /  R,(t,v)  co8  2irf0t  cos  2t/0  vdtdv 

Jo  Jo 


(E.l) 


where 


R.  (<>«)=/  /  ht(t-T)hi(v-'1)R,{T,'i)Re(r-'))dTd'i  {E.  2) 

J  —  oo  J  oo 

and  Rc(t  —  *■/)  is  defined  as  the  autocorrelation  function  of  the  spreading  code.  For 
DS-BPSK  spreading  modulation,  J2<.(6)  is  well  known,  namely 

Re(6)  =  (  1  ~  rf  ;  <  T< 

[0  ;  \6\  >  Te 

where  T~l  is  the  digital  rate  of  the  spreading  code  and  is  called  the  chip  rate. 
The  power  spectrum  of  the  spreading  code  Se(f )  is  the  Fourier  transform  of  Rc{6 ). 
namely  [Ref.  2]. 

Se(f)  =rcsinc2(/rc)  (£.3) 


where 


Sin  T x 


Figure  E.l  Autocorrelation  Function  and  Power  Spectral 
Density  of  DS-BPSK  Spreading  Code 


A  diagram  of  i2c(£)  and  Se(f)  is  shown  in  Fig.  E.l. 
A  new  signal  Pe  (t)  is  defined  as  follows: 


f  cos 2jt/0*  ;  0  <  t  <T> 
\  0  ;  otherwise 


So  that  equation  E.l  can  be  rewritten  as: 

a)  =  [  f  R,(t,v)Pe[t)Pe(v)dtdv 

J  —  00  j  —  oo 


(E.4) 


It  is  possible  to  express  R M  (t,  v)  as  the  double  inverse  Fourier  transform  of  Sz  (/,  v), 
that  is 

«,((,»)=/”  r  “'ifiv  (E. 5) 

J  —  oo  J  —  oo 

where 

Sz(f,v)=  r  r  R.{t,v)e-i2w{ft-’^dtdv  {E. 6) 

J  —  oo  J  —  oo 

Let  Pc(/)  be  the  Fourier  Transform  of  Pc(t),  so  that  Eq.  E.2  becomes 


°*j= r  r.  sz(f,u)  r  r  p^^^p^-^dtdvdfdu 

J  —  oo  J—  OO  j—  OO  «/ —  OO  Jj 

(/,  i/)Pc  (  — /)Pc 

Using  equation  E.2  to  evaluate  equation  E.6  results  in 

S*0»=  f°°  P,(r,7)Pe(r-7)  /"  /"  (t  -  r)e-'3"‘. 

J  —  oo  J  oo  «/  —  oo  J  —  oo 

ht(v  —  7)eja’,‘/*dtdvdrd7-  (£-8) 

=  Hr{f)Hj{-i/)  f  j  P,(r,7)i?e(r  -  7)drd7 


since 


and  Hj  (/0)  is  the  Fourier  TVansform  of  ht(t).  Therefore 


o>  =  r  r  HMH^p.i-np.w 

(£.10) 

A;(r.7)&(r  -  l)e-i2m'9 dr*i4fdv 

The  double  integration  involving  the  variables  r  and  7  in  equation  E.10  can  be  eval¬ 
uated  for  the  special  case  where  y(t)  is  the  nonstationary  random  process  described 
previously. 

In  order  to  evaluate  0* ,  the  double  integral  D(r,  7)  defined  as  must  first  be 
evaluated,  where 


D(r,~,)  =  J~  J~  J4(r,7)/2.(r,7)«-,a"Va'^drd7  (£.11) 

Recall  from  equation  2.20  that 


£,(»-, 7)  =  2cos2*/0(r -7)»j(m)  (£.12) 


where 

^(r*^)  =  j  *i(f-<*)*t(  1  -  a)q,(a)da  (£.13) 

The  double  Fourier  transform  pair  ij(r,  7)  and  H(u,w)  can  hence  be  expressed  as: 


£(«,«/)  =  J  J  ri(r,  ■y)e'  ,3wr*ej3ryw  drd-j 

so  that 

q(r, 7)  =  J  j  H(u,w)e)3rr,‘tj39ymdudu; 


(£.14) 


(£.15) 


va 


Therefore, 


D{rn)  =  J  /"  2[/~  /"  &{***)***”* *-*am‘,md*dw  • 

R,(t  -7)  coe2ir/o(r  —  ~j)t~ *7" 1 T e?lw‘'1  drdri 
~f  f  H(u’w)  [  f  2Re(T  ~  7)  cos  2jr/0(r  —  7)* 

J  —  OO  J  —  OO  •/  —  OO  */  —  OO 


-•)*  i2w(m  -  v)i 


drd'idudw 


=  r  r  ff(u,wj  r  [s.(/ -«+/„) +$.(/- 

J  —  OO  *  —  OO  J  —  OO 


■j3»(/ j3*(«  -Hi 


d'ldudw 


=  j  f  H(u,w)[Se[f  -  u  + f0)  +  Se(f -u  -  f0)]6(f  -  u  + w  -  i/)du 

J  —  OO  J  —  OO 

(£.16) 

Carrying  out  the  last  integration  with  respect  to  u;,  results  in 

£>(r,-r)=  I"  H(u,u  +  v  —  f)[Se(f  —  u  +  f0)  +  Se(f  —  u  —  f0)}  du  (£.17) 
J  -  OO 

Evaluation  of  £(u,w)  which  is  defined  in  Equation  E.14,  can  be  accomplished 
by  performing  a  double  Fourier  transformation  on  rj(r,^).  The  single  Fourier  trans¬ 
form  of  n(r,7)  is 

/  Tf(r,'i)e~ l2amr dr  =  ^  hL(^  -  a)  J  hL(r  -  a)e~}3wmT drq.(a)da  (£.18) 


Since 


[  hL  (r  -  a)e~i3muT  dr  =  («)«-'**' 

J  -  OO 


transforming  again  the  left  hand  side  of  equation  E.18,  results  in 

/_  r  rl(r'~l)e~’3*ure’3w'"'drd'1  =  /°°  J5M»)ff.(a)e- 

f  hL  (7  -  a)eJ,wwl  d~fda 
J  -  OO 


- j2wua 


H(u,w)  =  Hl(u)Hl(-w)  [  q.[a)e-}3 


wlu  —  w)a 


(£.19) 


'  y*,* V  v’V  vV'1?1 


However,  being  a  periodic  deterministic  signal,  qt  (a)  can  be  expressed  in  terms  of 
its  exponential  Fourier  series  expasion  as  follows 


I 


$ 


*vVl 

*jV? 


.VI 


«.(«)=  E  c-*"1 


so  that 


I  q.(a)c-’3w{u-m)ada  =  C,  f 

J  -OO  _  J  -OO 


(£.  20) 


Evaluation  of  the  C„  coefficients  for  a  particular  case  of  interest  is  worked  out  in 
Appendix  F.  Combining  Equations  E.19  and  E.20  yields 

OO 

H[u,w)  =  ^  CnHL{u)HL  {—w)6{nR,  -  u  +  u>)  (£.21) 


From  equation  E-21,  it  is  clear  that 


H(u,u  +  u-f)  =  HL{u)HL{f-u-u)  jP  C„*(n£,+i/-/)  (£.22) 

ns  -  oo 

Equation  E.17  can  hence  be  further  simplified  to  yield: 

f*  OO  OO 

D(r,'l)=  I  £  CnHL{u)HL{f-v-u)6{nR,+v-f)- 

J  —  OO _ _ 


[$>(/  +  /•-«)  +  Sc{f  -  fo  -  tt)]du 

=  E  C.H*R,  +  »  -  f)  r  BtWBtif  -  »  -  *)• 

f%m  —  oo  J-OO 

[S.if  +  fo  -«)  +  S.(/-/o  -u))du 


(£.23) 


Observing  that  in  the  integrand  of  equation  E.23,  u  is  the  variable  of  integration, 
we  can  define 

»L  (/•>')  =  f“  BL(u)BUf-u-u){S.{f  +  h-u)  +  S.{f-h-*)\du  (£.24) 

J  —  OO 


V,*- 


v«v 


Substituting  now  D(t,i)  in  equation  E.l,  with  the  use  of  equations  E.23  and  E.24 
yields: 


,r. 


*5“  T  r  B‘U)B,(-»)P.(-f)P.W)  t  C.((nR,+u-f)HHf,u)du df 

j  —  oo  J  —  oo 


(£.25) 


Integrating  Eq.  E.25  over  the  variable  u ,  we  have: 


t*  oo  oo 

«5  =  /  £  C.fl,  (f)S,  (nR,  -  f)P,  (-l)P.(S-  nR,  )H‘L(fJ-nR,)df 

j  “  OO  _ 


(£.26) 


From  equation  E.24,  H[  (/,  /  —  nJZ, )  can  be  seen  to  become: 


Hi{f,f-nR^)=  (  ££(u)£L(n£4-u)[5e(/+/o-tt)+Sc(/-/0-tt)]du  (£.27) 

J  —  OO 


so  that  finally 


<4  =  EC./"  (/-„«,). 

f  (u)Hl  (n£,  -  u)[S„(/  +  /0  -»)  +  £,(/-  /0  -  u)]du# 

*  -  oo 

(£.28) 
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THE  EXPONENTIAL  FOURIER  SERIES 
OF  q,{t)  AND  AVERAGE  POWER  OF  THE  PULSED  JAMMER. 


Given  q,  (f)  as  defined  in  Chapter  II,  part  C  and  re*illustrated  as  shown  in 
Figure  F.l  since  q,  (t)  is  a  periodic  signal,  it  is  possible  to  expand  q ,  (t)  in  terms  of 
an  exponential  Fourier  series  with  coefficients  C„ ,  where: 

C  >[  f'T'  AVa'"‘/T'dt+  /  '  {F.l) 

*•  L  Jo  J*r, 

Let  ti  —  t  —  ^  in  the  first  integral  and  let  t?  =  t  -  m  the  second  integral, 

so  that 

Cn  =  —  /  Aae,a**(*‘+#'^,/T*dtl  +  /  eae^s.Mi,+  -«-7 — i/r,^ 

(P.2) 

Using  properties  of  even  and  odd  functions  it  is  possible  to  simplify  Eq.  F.2  to 
yield 


=  ~ 

*1 


2A3 et2wH*l2  ("  cos?^^  +2CVa'",l">/a  /“  '*  *  cos  2-^-dt, 
Jo  T<  Jo  T, 


1 

9/iapi¥a«in2nxt1/rf 

»r, 

»  1 

•+■  — 

_j3  sin 2 nxta/Tf 

r. 

L  2mr/T,  j 

T 

0  M 

2nx/T, 

( i  -  #>r« 


sinn*p  « ,  .  sin  n*(l  -  p) 

=  AW" - -  +C3  1  -p  p 

nxp  nx(l  -  p) 


From  equation  F.3,  we  observe  that  if  n  =  0,  then 

C0  =  Aap  +  C3(l  -p) 

73 


(P.3) 


(P.4) 


Figure  F.l  Function  q,(t) 


In  the  cue  where  n  ^  0,  C„  can  be  reformulated  u 


C.  =  e3**' 


A2p 


sinnrp 

nwfi 


+  (-l)-C'(l-->) 


ainnir(l  -  p) 
nir(l  -  p)  J 


,n/0 


(F.5) 


However,  expanding  sinnx(l  -  p)  and  taking  into  account  that  sinn*  =  0  and 
coenjr  =  (—1)",  the  second  term  of  equation  F.5  becomes  -p3 C 3  ,  so  that 

we  finally  obtain 


CH  =  e*m00p3{A2  -  C 2) 


sinnirp 

nirp 


n  ^  0 


(F.6) 


Summarising  the  results  above,  we  have 


Cn  = 


I 


Aap  +  C"(l-p) 
p3(Aa  -  C* **J> 


if  n  =  0 
if  n  ?  0 


IF.  7) 


The  average  power  of  the  pulsed  noise  jammer  can  be  obtained  directly  from 
the  coefficients  C„ .  Alternatively,  recall  the  expression  for  A, (Ms)  from  equation 
C.I  and  let  tx  =  ta  =  t  so  that 


=  £,(*)  =  2 cos 2*/„(t  - 


(*  -  a)ht  (t  -  a)q,  (a)da 


or 


£,(*)=  2  f  h\  (1  -  a)q.  (a)da 
J  -  oo 
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(F.8) 


Assuming  specifically  that  hL  (f)  =  Bt~ B,U{t),Rt{t)  can  be  obtained  ss  follows, 
*,(<)  »  J  f~  {B,-‘<,-<U(t-a)]‘,.(a)da 

em  » 

=  2 B2  t~2Bm  V  C.e,a—“-',/r’dz  where  x  =  t  -  a 

Jo  _  _  __ 


=  2fla  £  C„e>a“* 


«/». 


e(3B  ♦ /3»«/T,  )» 


~{2B  +  j2*n/T') 

gjlmwt/T, 


=  2B1 


gjlWHt/t, 


2B1 .5.  C"  3fl  +  >2»»/T, 

C»  Y'  y.  F- . . 

2 B+  "2 B-  j2n*/T,  +  ^  C*  2fi  +  >2nir/Tt 


(F.9) 


=  2B7 

Since  the  Cn  coefficients  have  been  shown  to  be  of  the  form 


Cn=e>'"fn(A,C,,>) 


(F.  10) 


Equation  F.9  becomes 


=  2  B2 


=  2Ba 


gj(i*»»+3»««/rt ) 


+  (2nir/rf)3eJ“-"1  "'/»**• 


/-*  "  cos  ( *5^  +  n*y>  -  tan' 1  **-) 


JB  ■  v/(2B)»  +  (*»#/Tt)> 

(F.10) 

It  is  obvious  that  /2,(t)  is  periodic  in  f  with  period  Tn,  since  the  average  of  a  cosine 
signal  is  zero,  we  therefore  have 


*<*>-»(£) 


or 


Pj  =  R,(t)  =  B[A2p  +  C2{l-p)\ 


{F.  11) 
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APPENDIX  G 


VARIANCE  a2j  FOR  FAST  JAMMER 


In  the  cim  of  fast  jamming,  the  variance  a2  of  the  jammer  as  demonstrated 
in  Eq.  3.18  simplifies  to  the  following  expression 

=c„  r  ij,(/)i>ip.(/)i*  r  \HL(u)fs;(f  (g.i) 

J  -  m  J  —  9o 

where 

s;(x)  =  s.(x  +  /„)  +  5,(i-/0) 

In  order  to  simplify  the  notation,  denote  the  Fourier  Transform  and  the  inverse 
Fourier  Transform  operation  by  /’{(•)}  and  F~ 1  {(■)}  respectively  and  let  *  denote 
the  convolution  operation.  Define 


(G.2) 


where 

V(/)=  fm  \HL(u)\2S'e(f-u)du  (G. 3) 

J  -  00 

so  that  a2  can  be  reformulated  as  follows 

o’  =  Co  r  \Ht{f)\2\Pc(f)\2V{f)df 

J  —  90 

c.  r  \H,U)\’\P.(f) r  r 

J-m  (G.4) 

\H,(mP,(f)\*e  ‘'-''dt  df 

=c0  r  v(t)F-‘  {iff, (/)iai/»  (/)!>}  it 

J  -  90 
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Recall  from  Chap,  m,  Eq.  3.17  that: 


/>.(/)  =  G(1  -  /„)  +  G(/  +  /„) 

C(/)  =  |.mc  (/T,)e'',"T‘ 
where  for  reasonable  values  of  /0  and 

|-P«(/)r  =  Pc(f)PA-f)  »  I G{f  -  fo)\3  +  !C(/  4.  /o)!3  (G. 5) 


or  equivalently 

I /’.(/)!’  =  (y  V  lint’l/  -  /„)T.  +  (£)*  ainca(/  +  /.)T.  (C.6) 


The  inverse  Fourier  Transform  of  this  function  can  be  easily  obtained  from  tables 
(see  Ref.  9),  namely 

r‘{ie.(/)l*}- 7* (£)(*'’"•'  ic*) 


where 


,ti  <  r* 

i*i  > 


The  same  considerations  apply  for  jff,  (/))’  since  Ht  (/)  is  the  bandpass  filter  trans¬ 
fer  function  specified  in  terms  of  shifted  lowpass  equivalents.  Since  H,  (/)  is  signif¬ 
icant  for  /  =  ±  f0  over  a  restricted  frequency  band. 


where  the  lowp. 


equivalent  HL*(f)  is  given  by 

1 


*ts(/)  = 


*  *  ;2 »/ 


so  that 


"-(/)'*  2  (^»  -  (2w/)>  ) 
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((7  81 


((?  9^ 


it;,  to) 


Consequently, 


F~l  {| Ht(f)\a}  =  [e'3"**  +e-'3"“] 


(C.ll) 


so  that  evaluation  of  F~l  {j/f/(/)|3|Pc(/)|3}  can  be  accomplished  by  using  the 
convolution  theorem  [Ref.  0],  that  m 

r-'  {\B,W\P.(t)P)  ~  F-'  {'B, (/)!*}.  F-‘  {/>(/)!’} 


-/:‘f‘(s)c . 

=  ^  j~  «•"-■  a(^)co.J»/.(I -t)it 

-fiy . ‘few 


g-  ^  g- H«- » 


1  2  w  f  t 


G  .12' 


where 


f(l.r)  =  coe2v/0r  coa2ir/»(t  -  r) 

!  131 

=  -  xos2*/0t  -*■  cas2v/0*coe4«70r  tin 2w /0 fsintvA,  r 

According  to  the  value  of  t  -  r  ,  Eq.  G.12  can  be  further  expressed  as  follows 

p  1  { #r,(/v3  p.(/v3}  =  /,-/,  k;  h 

where: 


/ 

‘  •  j 

r.*  “■ 

G  IS 

/ 

*  *  J 

r-  "x(k)*,-,ur 

16 

We  can  now  proceed  to  the  specification  of  r|f’  which  a  obtained  from 
*(l)  =  F -‘ly  ff ■  (w)  ‘  5*{/  elds  | 


r  '  {  hl  if)  * }  r  1  is  i/u 

a#  9  ‘  (  p  )  c°e  2v/.f 


;  1 ' 
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since 


F~  l{|#t  (/)|a}  =  Be”*"1  eoa2»/0t 


Substituting  Eq.  G  .4  with  Eq.  G.14  and  G.17  yields 


(G.  18) 


=C0B  J~  e-*'*,A^)coe2ir/0<(/l  -  I7)dt 

The  limit  of  integration  can  be  reduced  to  the  interval  -T..T,  since  the  function 
A*/0t(/i  A )dt  The  limit  of  integration  can  be  reduced  to  the  interval  -Te,Te 
since  the  function  A  ^  ^  is  sero  except  for  <  <  Tt .  hence 


-  BC  J  e  ■*  *  A  (  ^  )  cos  2» /..<!/  ~I:dt 

According  to  the  value  of  t,  two  caaes  can  be  considered  u  follows 
Case  (1)  T*  v  t  v  o. 

In  this  case,  we  have 


(7  19 


/  -  /i  = 


f  '  “*  1  *  jr  )  f {t.r)dr  +  J  t  * "  ’  (  1  T  )  g\t.  - 


idr 


•r**  ■('  f.)'"’ 


Idr 


(7  20 


Since  Eq  G  13  shows  that  fjl.r)  «  «  function  of  cos  4v /, »  and  sin  4*  *  it  can  be 

seen  that  the  integration  of  functions  of  the  form  f** '  cos  4v  /.« r  and  e“* '  nn  4»  ' 
as  well  as  forms  If"'  cos4»/.'  and  If*'  sin4v/,t  will  vield  terms  involving  the 
actors  m  -  I6v:  cl  and  *n*  *  16v;  1  :n  the  denominator  >>t  the  *es  i  ’  vvr  .,r 


the  numerators  remain  bounded  The  numerators  are  proportion*,  *o  m 
this  case  equals  or  f?  Since  in  a  practical  system. 


ich  T 


fc  and  ( 


AT 


/I 


the  contribution  of  the  corresponding  terms  are  therefore  not  significant,  so  that, 
Eq  G.20  can  be  simplified  as  follows: 


h  +  I*  «  ^  •  ^cos2ir/0t  je~6‘  J  ebT  ^1  +  dr+ 


(G.21) 


Eq.  G.21  can  be  evaluated  by  using  standard  mathematics  tables  [Ref.  10]  and 


therefore 


,  r  „  ,  1  t  cbt  e~fcT‘  ,  bt  bt 

h  +  h  «  ~ cos 2jt/0 t  -  -  — —  — —  +  — — —  (e~6  +  e 

8  b  bTh  b?Th  2 PTb  v 


Case  (2)  0  >  t  >  -Tb, 


(G.  22) 


In  this  case  we  have 


= t  m  (■ + y  j(,’r) + /  (* + s)  a{t'r)dT 

♦r- . 

(G.23) 

By  considerations  similar  to  those  presented  in  Case  (l),  Eq.  G.23  can  be  approx¬ 
imated  by 

T '  i /r> **r  (l  +  £)  *  + /  *■*'  (» +  £)  * 


(G.24) 


so  that 


,  r  W*  „  ,  Tl  *  C6‘  6“^ 

/.  C0.2./.I  i  +  --  —  +  — 


(e~bt+ebt)  (G.25) 


>  *  v 


.lt  *rT|.;T| 


We  can  now  proceed  to  the  evaluation  of  a2  expressed  by  Eq.  G.19  according 
to  the  range  of  t.  Since  Te  <  Tb, 

BCobTb  Bt(  t\  \1  t  S*  (  bt  bt 

^  “  8  j_Te  v1+  tJ  2  fot[b  +  bTb  (PTb  +  2lPTb{  + 

BC0bTb  [T°  (  _  £\  3  fl  _  _L  _  £_  (  h  bt 

S  J0  V  te)  fo[b  bTb  PTb  2b>Tb\ 

{G.  26) 

and  using  cos2  2it f0t  =  j(cos4 irf0t  +  l)  the  reasoning  that  gave  rise  to  the  simpli¬ 
fications  in  evaluating  Eq.  G.20  can  be  used  here  to  neglect  the  terms  involving 
(m2  +  16k2/2)  or  (m2  +  16ir2/^)2  in  the  denominator.  Thus  Eq.  G.26  can  be 
further  simplified  as  follows: 

2  BCobn  r  r  .a*  a  .  *  \  ri ,  *  *6t  ,  ( ~~bt ,  ,Mn1  ^ 

7  16  1  J_T  l1+  Tj  b  +  bTb  fPTb  2b2Tb  (  ' 


(G.27) 


t 

ebt 

1 

e~hT * 

bTb 

lPTb 

ebt 

e-bT> 

I 

(«-  ■ 

tPTb 

2b2  Tb 

and  using  standard  mathematics  tables  (Ref.  10)  results  in 


ai  = 


=  BoCoW;  flTj.  __1_/  _  Br  \1  J_r_2_/  B 

'  8  \b[B  B2TC  V  ;J  +  6rfc[B3rcV 

1  (1  I  1  [  C~Br0g*r‘  I  g~BT°g'6rc 

B2\  J\  PTb[2[B -b)2Te  2{B  +  b)2Tc 

B  B2+b2  I  e-6r*  f  2 B  2  (B2  +  b2) 


B2  —  IP  {B2-b2)2Tc\  lPTb  [B2  -  IP  ( B 

g-(B-b)Tc  e-(B  +  b)T.  -M 

+  (b-B)2Te  +  (B  +  6)2rJ  J 
which  by  regrouping  terms,  obtains  the  final  form 


(G.28) 


B(B2  +  63) 


bTb(B2-(P)  bTbTc(B 


Bc-BT‘e~bT' 

Be-BT'cbT<  ' 

) 

2bThTr  (B  4-  b) 

2  2  bThTc{B-b) 

2 

1 

STATISTICS  OF  Y,  AND  St 


th 


We  first  recall  the  definitions  of  Nih  and  Y}  from  Chapter  IV  as  follows 

rn 


Yj  =  /  y(t)e(t)  *  ht  (t)  co»2xf0tdt 
Jo 

N*k  =  f  Nth[t)c(t)  •  ht(t)  cot  2* fatdt 
Jo 


(H.l) 

{H.2) 


Observe  that  Yj  is  a  zero  mean  random  variable  since  the  jammer  and  the  spreading 
code  c(t)  are  assumed  uncorrelated  and  furthermore,  e(t)  has  zero  mean.  The 
variance  a*  of  Yj  has  been  developed  in  Chap.  3  and  can  be  expressed  as  (see  Eq. 
3.19) 


2  _  CoT> 

Oj  -  — 


BT. 


8 


+  (1  _  t~BT')( _ - _ —  ^  -  - till 

\  ’{B'TtT'  BTj  BT> 

"[i m--' 


Be-*T't-hT‘ 

+  _ _ _ rrr  + 


ft3)  6Tbrc(B3  -vy 

Bc-BT'ciT • 


i,H.  3) 


} 


2bT„Tc (B  +  6) 3  2bThTc (B  -  6) 3 
Replacing  now  B  and  b  in  Eq.  H.3  by  2 *RC  and  2-kR^  respectively  thereby  setting 


the  3  dB  cutoff  of  the  filters  appropriately,  with 


_  4  Re 
G--% 


defined  as  the  processing  gain,  a  final  expression  for  <r3  is  obtained  as  follows 

=  ^J(G,)  (ff.4) 

where 

1  ,  ,  (1-«-]')(1-tG.)  1 

2  xG, 


,(r  ,  _  ,  ,  (1 -«-»■)(!-*(?.)  _  !  +  «-»-  -  1)  f  1 

'  c'  n—a/~ i  o_n  r»_  j 


2ir3G 

e 


e 

-2w 


r 


2n 

n 


(H.5) 


Gl  +  1 


_  J_JL 

e  a* 


li 


The  determination  of  the  variance  °'tK  of  the  random  variable  Nth  can  be  ac¬ 
complished  in  a  manner  very  similar  to  that  used  in  the  determination  of  o* .  Ob¬ 
serving  however  that  replacing  in  the  jammer  model  the  deterministic  priodic  signal 
q(t)  with  a  constant  unit  amplitude  function,  namely  q(t)  =  1,  setting  Rw 
equal  to  ^6 (see  Eqs.  2.16  and  2.17)  results  in  a  simple  AWGN  noise  model 
having  PSD  level  N0/ 2.  Thus  all  evaluations  involving  the  jammer  can  now  be  used 
with  the  stated  modifications  in  order  to  characterize  the  noise  contributions  to 
the  performance  of  the  system.  Therefore,  (see  Eq.  E.13) 

*l(r»7)  =  /  (r  -  a)hL  (7  -  a)da  {H. 5) 

J  —  OO 

and  with 

hL(t)  =  Be~BtU{t) 


then 


r?(r^)  =  J' 


Mln(r.ir) 


e2  B  “  da. 


~  2 

Thus,  the  auto  correlation  function  R„  ( t ,  r)  of  the  noise  Nth  (t)  is 


(H.6) 


R.(t,r)  =  yBe'®1"11  cos2jr/0(r  -  -y) 


(H.7) 


From  the  methodology  of  Chap.  3,  Eq.  E.14  becomes 


H(  u,tv)  =  f  j 

J  —  OO  *  *  OO 


(H.  8) 


Evaluation  of  Eq.  H.8  can  be  shown  to  yield 


H{u,  w)  =  Nt 


-N*i 


+ 


2(5  —  j2iru  B  + 


(H.9) 
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Therefore,  the  variance  afh  of  Nth  (see  Eq.  E.25) 

f  H{u,u  +  u  -  f)[Se{f  -  u  +  fa)  +  Se(f  -  u  +  f0)]dudfdu 

J  —  oo 

(jy.io) 


where 


g(«,»+ ^  -  /)  -  w.f  {  a*  -hL)*  }{(“-(u  +  ‘<- 

=jv»f{^+(Lp}<(/-‘/) 


(g.n) 


Letting  /  =  u  —  fQ  or  /  =  u  +  fa  as  appropriate  in  the  two  integrals,  obtain 

<= r  jv.ifl'I(/+/.)n(/+/.)ij^(/)<i/+r  H'\H<v-f')p.u-h)?z{M 

J  —  OO  J  —  oo 

(H.  12) 


where 


/oo  p  OD 

..TsrnssF**'-)**  <*13> 


Z(/)  can  be  evaluated  further  as  follows 
B  2  B  _  B 


ZW=2W( 2^  **•<'>  =2F{e'*,,|  *‘<‘>} 


_b  r 

■  2  y_0 

-*ir 


e-fl(l,iZa(t)e"y3'/‘^ 


(F.14) 


e”fl‘  cos  2irftdt 


-r«-^ 


cos  2irftdt 


The  two  integrals  in  Eq.  H.14  can  be  evaluated  separately  and  with  STC  =  27r,  we 


*(/)-sr(/> 


where 


z'(/)  =  [1  + Jrj,,,  [(«•”  -  i)  ((/rc)a  -  i)  -  2 /r.  «n2x/r.) 

(g.15) 


From  Figure  H.l,  it  can  be  obeerved  that  if  <<  R,t  then  Z'(f)  w  1  for  /  within 
the  interval  In  fact,  Z(f)  ia  significant  only  for  |/|  <<  H*.  Hence,  Eq. 

H.12  can  be  simplified  ss  follows: 

<«•£[/*'  \H,U  +  f.)P.V  +  h)l'<V  +  /*'  l«, (f-f.)P.lf 
=  v/*  l*M/ -/.)*(/ 

*  .'-it. 

(//.16) 

From  Eq.  G.6,  we  have 

ift(/  -  /.)|*  *  |G(/  -  /  -  */.)!’  +  |G(/)|J 

where 

<?(/)  =  (j)  >inc2 fT„.  (H.17) 

The  term  \G(f  —  2/0)|  has  significant  components  outside  the  integration  region 
[— fZfc.fZfc].  Furthermore  from  Eq.  G.8 

\B,(f  -  /.)!’  *»  |S..(/ -2/.)|a  +  l^.(-/)r 

and  again,  the  term  (/  —  2/0)|  can  be  dropped  from  further  consideration  due 
to  reasons  similar  to  those  discussed  above  in  the  context  of  \G(f  —  2fa)\.  Thus 

<  »  (D’sincVr.d/ 
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and  since  \Htt  (/))*  has  3  dB  cutoff  at  /  =  ±R+,  the  integral  itself  can  be  approx¬ 


imated  by 


Finally 


■inca(/Tk)tf  «  *  =  0.903 


'tk 


S3 


Mt 

2* 


k 


{HAS) 
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